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INTRODUCTION 


The  purpose  of  this  study  is  to  develop  analysis  techniques  for  blast 
loaded  parabolic  antennas  usinq  existing  structural  codes,  such  as  NASTRAN. 
The  analysis  techniques  are  for  application  to  a UB  Army  Satellite  Com- 
munications Agency  (SATCOM)  Antenna  to  determine  its  behavior  when  sub- 
jected to  air  blasts.  The  antenna  consists  of  a five  niece  monocoque 
reflector,  lightweight  reflector  backup  castings  for  the  gimbal  and 
actuator  attachment  points,  a cross  elevation  axis  of  rotation  with  linear 
actuator  drive,  a tubular  elevation/cross  elevation  gimbal  assembly,  an 
elevation  axis  of  rotation  with  linear  actuator  drive,  and  a double  tripod 
supporting  the  elevation  axis.  The  problem  is  simplified  and  streamlined 
wherever  possible  to  avoid  excessive  computer  running  costs.  Computer 
programs  are  described  which  reduce  the  blast  parameters  and  antennas 
geometry  to  concentrated  loads  at  the  NASTRAN  grid  points.  Provisions  were 
made  to  account  for  nonlinearities  related  to  resistance  to  certain  grid 
point  displacements. 

The  objective  of  this  study  was  to  determine  gauge  locations  and  a 
nondestructive  load  range  for  the  SATCOM  antenna  when  subjected  to  frontal 
shock  tube  blasts.  The  conclusions  stated  here  will  be  based  on  the 
geometry  of  the  antenna,  a NASTRAN  static  analysis  of  a model  of  the 
re<^lector,  and  a report  by  Tasker  Industries  on  finite  element  studies  for 
nuclear  hardness  of:  a)  the  antenna  on  its  shelter,  and  b)  the  antenna 

center  section  dismounted  from  the  shelter. 

DISCUSSION  - Analysis  of  a Blast  Loaded  Parabolic  Antenna 

Although  the  antenna  system  is  not  structurally  symmetric,  the 
relatively  high  stiffness  of  the  support  in  comparison  to  the  reflector 
makes  it  possible  to  conservatively  analyze  the  reflector  as  if  it  were 
fixed  at  its  connection  points  to  the  supports  and,  therefore,  symmetric 
about  a horizontal  axis.  This  assumption  is  expected  to  reduce  the 
computer  running  time  considerably,  but  has  the  disadvantage  of  restricting 
the  model  to  the  study  of  symmetric  loads.  It  should  be  noted  that  avail- 
able analyses  indicate  that  a head-on  (symmetric)  blast  loading  will  be  the 
most  critical.  The  NASTRAN  model  for  the  reflector  is  shown  in  Figure  1 
through  3.  Membrane  plates  were  used  to  model  the  skin  and  ribs  of  the 
petals  and  center  section,  except  at  connection  points,  where  membrane- 
bending Plates  were  used.  An  alternate  model  usinq  all  membrane-bending 
plates  makes  it  possible  to  study  the  flexual  stresses  to  some  extent. 
Elastic  springs  were  used  to  model  the  rim  latches,  the  tx'ta  1 -p«'ta  1 latches, 
and  the  petal-center  section  latches.  Nonl inear it ies  exist  where  the 
petals  join  together  and  where  the  petals  join  with  the  center  section, 
klhen  the  joining  parts  are  pressed  together,  a resisting  force  is 
developed;  when  these  points  movi'  aoart,  no  resisting  force  develops.  An 
apF>roach  wa-s  devi;  d using  nonlinear  springs  to  account  for  this  behavior 
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in  the  transient  analysis  only. 

The  mechanics  of  this  approach  is  explained  in  Appendix  A.  It  is 
recommended  that  the  nonlinear  springs  be  omitted  unless  there  is  an 
indication  that  the  nonlinear  effect  will  be  siqnificant.  The  support 
points  and  most  of  the  latches  on  the  reflector  are  cast  aluminum.  These 
were  not  modelled  precisely  because  accurate  dimensions  were  not  available 
and  existing  analyses  indicated  a precise  model  was  not  necessary. 

However,  changes  were  made  in  the  reflector  model  in  the  vicinity  of 
latches  and  supports  to  approximate  their  influence  on  the  rest  of  the 
structure.  Details  of  the  reflector  model  in  the  vicinity  of  the  petal- 
petal  latches  and  in  the  vicinitv  of  the  petal-center  latches  are  shown  in 
Figure  4.  Details  of  the  reflector  supports  for  the  symmetric  model  are 
shown  in  Figure  5.  The  cross  elevation  actuator  is  included  in  the  model 
and  is  assumed  to  be  pinned  at  its  support  areas,  while  the  gimbal  support 
is  constrained  for  all  motion  except  rotation  about  a vertical  axis.  It 
was  convenient  in  most  cases  to  define  grid  point  loc=>*-ions  and  movements 
in'  terms  of  a cylindrical  coordinate  system  with  origin  at  the  center  of 
the  back  panel.  Plates  were  generally  defined  so  that  outputted  stress 
directions  would  be  in  the  meridional-tangential  directions  for  elements  on 
the  two  faces,  in  the  meridional-axial  directions  for  ribs,  and  in  the  axial- 
tangential  directions  for  edge  elements.  If  an  analysis  is  desired  for  the 
center  section  alone,  all  grid  points,  elements -and  load  cards  related  to 
the  Petals  and  all  elastic  springs  must  be  removed.  Also,  the  case  control 
deck  and  the  bulk  data  deck  must  be  purged  of  all  nonlinear  loads  and 
multipoint  constraints  which  refer  to  the  petals. 

NASTRAN  has  two  possible  formats  for  describing  transient  loads.  One 
is  defined  by  the  conditions  P(t)=0  (f<0  or  T>T2-T1) 

P(t)=AtS  e ct  j-os  (2-fTFT+G)  (o<  t JET2-T1) 

where 

P(t)=the  time  dependent  load  at  a point  for  a degree  of  freedom 
t=t-Tl-  TT 
t=elapsed  time 
*f=time  delay 

B,C,F ,G,T2,T1  are  constants  defined  bv  the  user. 

The  other  NASTRAN  loading  format  is  of  the  form 

P(t)=A  F(t-T)  where  P(t),t,  and  T are  as  defined  previously 
F(t-T)  is  a tabulated  function  of  time  that  is  linearly  interpolated 
between  entries  and  A is  a tabulated  coefficient  which  mav  be  different  for 
each  loaded  degree  of  freedom. 

The  formula  approach  was  considered  first  because  of  its  apparant 
simplicity,  but  was  rejected  because  of  inconveniences  it  engenders  when  a 
chance  in  the  pressure-time  relationship  is  desired.  The  tabular  approach 


5 


is  applied  to  the  antenna  as  follows:  a pre-processor  program  generates 

the  tabular  information  defining  the  pressure  on  the  surface.  Only  a few 
of  these  points  will  generally  be  needed  to  define  the  time-pressure  curve 
accurately  in  the  NASTRAN  table.  This  orogram  is  described  in  Appendix  B. 
The  tabulated  coefficient.  A,  accounts  for  the  area  associated  with  a 
uniform  pressure  load  on  the  reflector  surface,  and  for  the  geometry  of  that 
surface.  A FORTRAN  IV  computer  program  was  written  which  reads  the  cards 
defining  the  front  or  back  surface  of  the  reflector  and  punches  the  NASTRAN 
DAREA  cards  which  define  the  constant.  A,  for  each  loaded  degree  of  freedom. 
Use  of  this  program  is  detailed  in  Aopendix  C.  The  approach  used  in 
splitting  up  the  area  is  the  same  as  that  used  by  NASTRAN  for  static 
pressure  loads.  For  triangles,  one  third  of  the  total  area  is  associated 
with  each  corner.  Each  diagonal  of  a quadrangle  is  used  to  split  the 
quadrangle  into  two  triangles.  The  triangle  areas  are  distributed  as 
described  above  and  the  area  associated  with  a corner  of  the  quadrangle  is 
the  sum  of  the  triangle  areas  associated  with  that  corner  divided  by  six. 

For  the  transient  analysis  of  the  reflector,  either  the  NASTRAN  direct 
or  the  NASTRAN  modal  approach  could  be  used.  The  direct  approach  is 
recommended  here,  primarily  for  its  ability  to  handle  nonlinearities.  It 
is  suggested  that  a static  analysis  and  a normal  mode  analysis  be  executed 
before  the  more  expensive  dynamic  analysis  to  check  for  weaknesses  in  the 
model.  These  preliminary  analyses  may  also  provide  insight  on  use  of  the 
Guyen  Reduction,  which  can  reduce  processor  time  by  decreasing  the  number 
of  points  given  inertia  properties.  The  time  increment  size  can  be 
roughly  determined  by  the  equation; 

At=  1 

’Of  where  A t=  time  increment  (secs) 

f=  hiqhest  frequency  expected  to  influence 
solution  significantly  (cycles/sec) 

The  Tasker  report  describes  a modal  transient  analysis,  using  the 
STARDYNE  finite  element  program.  Computer  processor  time  was  in  excess  of 
three  hours  for  the  analysis  of  the  antenna  on  the  shelter,  so  several 
measures  were  taken  to  increase  the  efficiency  of  the  analysis.  One  such 
measure  involved  minimizina  the  number  of  monitored  points,  to  the  extent 
that  only  latches  and  two  rib  plates  were  monitored  on  the  reflector.  Pre- 
sumably these  points  were  found  to  be  critical  in  a less  expensive  stafic 
analysis.  The  transient  analysis  considered  only  elastic  behavior  of  the 
structure;  since  time  dependent  effects  on  material  properties  were  either 
unknown  or  not  accounted  for,  the  Tasker  report  defined  failure  in  terms 
of  the  ultimate  strength,  i.e.,  the  stress  at  which  separation  of  molecules 
occurs  in  a static  situation.  A failure  factor,  determined  by  dividing  the 
maximum  stress  on  a member  by  the  ultimate  stress,  was  used  to  indicate  the 
severity  of  stress.  'Five  monitored  antenna  elements  "failed"  in  this  test. 
These  elements  and  their  failure  factors  are  shown  in  Table  1. 
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TABLE  1 

Suiranarv  of  Antenna  Failures 


MEMBER 


FAILURE  FACTOR 


Latch  Casting  6.1 
X-Elevation  Actuator  3.6 
Elevation  Actuator  7.2 
Gimbal  Tube  1.6 
Rear  Leg  1 . 4 


The  high  stress  level  in  the  latch  casting  is  indicative  of  the 
stresses  in  the  unmonitored  portion  of  the  reflector,  where  it  is  likely 
that  failure  occurred.  The  possibility  of  a buckling  failure  is  high,  but 
is  not  reflected  in  Table  1.  For  this  analysis,  the  shelter  was  available 
to  absorb  considerable  energy  and,  thus,  reduce  the  blast  effect  on  the 
antenna.  The  antenna  will  be  ground  mounted  for  the  shock  tube  test. 


In  the  second  analysis,  the  center  section  was  ground  mounted  and 
subjected  to  the  same  blast.  The  smaller  model  allowed  for  more  points  on 
the  antenna  to  be  monitored.  The  monitored  elements  which  "failed"  and 
their  failure  factors  are  given  in  Table  2. 


TABLE  2 

Summary  of  Center  Section  Failures 


MEMBER 


FAILURE  FACTOR 


Actuator  trailing  arm 
Rear  leg 

Reflector  skin  (worst  case) 


2.6 

1.2 

6.4 


The  report  indicates  that  every  reflector  plate  element  monitored 
experienced  stresses  exceeding  the  allowable  buckling  load  in  the  center 
section  test. 

A study  of  the  geometry  and  the  static  NASTRAN' solution,  as  well  as  the 
Tasker  report,  indicated  that  the  upper  portion  of  the  reflector  away  from 
the  cross  elevation  support  is  where  the  largest  stresses  are  likely  to 
occur.  This  is  because  the  stiffness  is  lower  and  because  fewer  paths  are 
available  for  removal  of  energv  from  this  area. 

CONCLUSIONS 

Analysis  of  a Blast  Loaded  Parabolic  Antenna 


A physical  model  and  blast  analysis  technioues  have  been  developed  for 
a solid  surface  pctrabolic  antenna,  using  the  NASTRAN  structural  analysis 
program.  Pre-processer  programs  were  developed  to  define  the  blast  pressure 
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load  in  terms  acceptable  to  NASTRAN.  A method  was  developed  for  dealinq 
with  nonlinear  deformations  cf  the  reflector.  It  is  expected  that  the 
model  will  give  reasonable  results  for  the  overall  behavior  of  the  reflector. 
If  information  concerning  localized  effects  is  desired,  a separate  detailed 
local  analysis,  based  on  the  solution  for  the  total  reflector,  can  be  made. 
This  separate  analysis  would  be  approximate,  but  should  be  sufficiently 
accurate  to  describe  the  localized  behavior.  Application  of  the  Guyen 
Reduction  is  recommended  for  subsequent  runs  once  the  general  dynamic 
behavior  of  the  antenna  has  been  established. 

Instrumentation  for  SATCOM  Antenna  Shock  Tube  Test 

The  Tasker  analysis  were  done  for  a specified  nuclear  blast.  Failure 
factors  are  given  for  some  of  the  members  which  failed,  presumably  the  most 
serious  failures.  A gross  estimate  of  the  maximum  nondestructive  over- 
pressure could  be  had  bv  dividing  the  Tasker  overpressure  by  the  largest 
failure  factor.  This  method  ignores  plasticity  and  buckling,  but  these 
phenomena  may  have  reduced  significance  due  to  the  high  strain  rate  and  the 
short  load  duration.  This  approach  does  not  account  for  the  energy 
absorbing  effect  of  the  shelter  which  was  included  in  the  Tasker  analysis. 

It  is  conservative  in  assuming  a linear  relationship  between  stress  emd 
overpressure  rather  than  stress  and  impulse,  since  the  duration  of  the 
nuclear  blast  is  probably  much  higher  than  that  of  the  shock  tube  blast. 

At  any  rate,  the  estimated  overpressure  should  be  used  as  an  indicator 
only. 

Instrumentation  was  concentrated  on  the  top  of  the  antenna  away  from 
the  cross  elevation  support.  These  locations  are  tabulated  in  Table  3. 
Location  codes  1,  2,  3,  4,  5,  7,  9,  12  and  16  indicate  points  where  high 
stresses  are  exoected.  The  other  p>oints  are  monitored  primarily  to  gain 
information  about  the  behavior  of  the  antenna. 
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TA3LE  3.  SENSOR  SELECTIONS  AND  LOCATIONS  ON  ANTENNA  FOR  SHOCK  TUBE  TEST 


Test 

W/Petals 

Location 

Code 

1 


Test 

W/0  Petals 

Location 

Code 


Sensor  Type 
Linear  Strain 


Channels 

Required Location 

1 Upper  center  latch,  left 

petal  (on  single  piece  center 
selection) 


2 


Rosette  Strain 


3 On  front  of  petal,  vicinity 

of  #1 


3 


4 . 


5 


Linear  Strain 


Linear  Strain 


Linear  Strain 


1 On  shaft  of  elevation  actua- 

tor (two  gages  in  series,  one 
on  each  side)  SEE  NOTE  1 

1 Rear  j.eg,  left  side  (two 

gages  in  series,  one  on  each 
side)  SEE  NOTE  1 

1 Gimbal  tube,  upper  left  side 

.(two  gages  in  series,  one  on 
each  side)  SEE  NOTE  1 


6 


7 7 


Linear  Strain  1 Left  petal  rim,  near  rim 

latch  connecting  to  upper 
petal 

Rosette  Strain  3 Between  ribs  in  center 

section  at  left  of  top  center- 
line 


8 


9 


Accelerometer 


Linear  Strain 


1 Upper  petal  rim,  near  rim 

latch  connecting  to  left 
petal 

1 On  shaft  of  cross-elevation 

actuator  SEE  NOTE  1 


10 


Rosette  Strain 


3 On  rear  of  petal,  vicinity 

of  #1 


11 


Rosette  Strain  3 On  front  of  left  petal, 

vicinity  of  upper  to  left 
(petal  to  petal)  latch 


j 

1 


a 


4 


9 


Test  Test 

W/Petals  W/0  Petals 
Location  Location 

Code Code 

12 


13  13  Linear  Strain  1 On  front  left  side  of  center 

section  over  and  perpend- 
icular to  rib 

14  Accelerometers  2 On  front  of  both  sides  of 

upper  petal  to  left  petal 
joint , away  from  petal  to 
petal  latch 

15  Accelerometers  2 Both  sides  of  petal  to  center 

section  joint  on  upper  part  of 
left  petal 

16  Rosette  Strain  3 On  front  of  center  section. 

Vicinity  of  #1 

17  Rosette  Strain  3 On  rear  of  center  section. 

Vicinity  of  #1 

18  Accelerometer  1 On  front  edge  of  center 

section,  near  petal  to  petal 
joint 

19  Rosette  Strain  3 Physically  on  edge  of  center 

section,  between  petal  to 
center  section  latches 

NOTE  1:  Location  #3,  4,  5 & 9 are  to  be  installed  on  structural  and  mech- 

anical members  which  support  and  position  this  antenna.  Directions 
as  to  their  exact  locations  will  be  supplied  at  the  test  site 


Sensor  Type 
Rosette  Strain 


Channels 

Required 

3 


Location 

On  rib  on  upper  side  of  left 
petal 
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APPENDIX  A 


NASTRAN  Model  of  Nonlinear  Behavior  of  Surfaces  in  Contact 

The  following  model  of  the  nonlinear  behavior  of  surfaces  in  contact 
for  the  reflector  should  be  used  only  when  this  behavior  will  have  a 
significant  effect  on  the  problem.  For  the  present  model,  a minimum  of 
thirty  eight  (38)  transfer  functions  and  seventy  six  (76)  nonlinear  loads 
must  be  developed  if  the  nonlinearity  is  to  be  handled  automatically  by 
NASTRAN.  The  procedure  is  as  follows: 

1)  Define  extra  points  for  each  pair  of  points  in  contact,  using  EPOTNT 
cards . 

2)  Relate  the  extra  points  to  the  points  in  contact  through  transfer 
functions,  using  TF  cards.  The  transfer  function  is: 

(1>  + (1)^1  + (-1)^2  = 0 where  Mu  is  the  extra  point.  For 

petal  to  petal  contact  surfaces  is  the  tangential  deflection  of  the 

first  of  the  contact  points  encountered  when  traversing  the  reflector  in  the 
tangential  direction  andZ^2  is  the  tangential  deflection  of  the  second  of 
the  contact  points.  For  petal  to  center  section  contact  surfaces,  ^2^  is 
the  radial  deflection  of  the  first  of  the  contact  points  encountered  when 
traversing  the  reflector  in  the  radial  direction  and  <^2  i®  ih®  radial 
deflection  of  the  second  of  the  contact  points.  The  presence  of  a transfer 
function  must  be  indicated  in  the  case  control  deck. 

3)  Apply  the  negative  power  function  nonlinear  load  to  both  points  in 

contact  using  N0LIN4  cards.  The  two  loads  are  defined,  respectively,  by 
the  equations:  ?2(t)  = -S^  (qUo  (t))l  , (t)=  -S2(qi<D  (t)<0 

P^  (t)  = P2(t)=0 

where  5^=  const  - arbitrarily  taken  to  be  the  same  as  used  in  the  elastic 
springs  and  82=  “Sj^.  The  presence  of  a nonlinear  load  system  must  be 
indicated  in  the  case  control  deck. 


APPENDIX  B 


COMPUTER  PROGRAM  FOR  PRESSURE 
TIME  HISTORY  OF  BLAST  LOAD  ON  THE  REFLECTOR  SURFACES 

For  load  definition  purposes,  the  reflector  can  be  compared  to  a closed 
rectangular  box,  a shape  which  has  been  subjected  to  extensive  study.  The 
empirical  formulas  presented  by  Gladstone^  for  the  rectangular  box  were  used 
in  a short  computer  program^  written  for  use  on  the  Hewlett-Packard  9830 
Calculator  to  provide  the  average  pressure-time  history  of  the  blast  load  on 
each  face  of  the  reflector.  Data  is  inputted  by  the  following  statement: 

9000  DATA  A,  B,  C,  D,  E,  F where 

A - Surface  indicator  (1  = front  2 = side  3 = back) 

B = Positive  phase  duration  (seconds) 

C = Time  increment  for  calculating  pressure  (seconds) 

D = Peak  overpressure  (psi) 

E = Length  of  Reflector  (ft) 

F = Heighth  of  "rectangle"  or  half  of  width  of  "rectangle",  whichever  is 
smaller  (ft)  (For  reflector,  F = radius,  in  ft) 

The  output  for  this  program  is  a printout  of  the  pressure  for  each  time 
step  and  a plot  of  pressure  vs.  time.  The  plot  is  used  to  determine  how 
many  of  the  printout  values  are  sufficient  to  accurately  define  the  load 
history.  Appropriate  data  can  then  be  selected  from  the  printout  to 
describe  the  load  on  NASTRAN  TABLED  cards. 

Formulas  for  load-history  descriptions  from  other  sources  can  be 
substituted  for  the  Gladstone  formulas  in  the  program. 


iThe  Effect  of  Nuclear  Weapons,  Samuel  Gladstone,  US  Atomic  Energy 
Commision,  Washington,  DC,  April  1962 

^Computer  program  and  plot  appear  on  following  pages 


10  READ  A,X1,D1,P0,L,S 
20  R=2*P0* (102.9+4*P0)/(102.9+P0) 

30  U=1117* (l+6*P0/7/14.7)  (1/2) 

40  SCALE  0,Xl*7/6,0,0,R*6/5 
50  OFFSET  Xl/8.5,R/5 
60  XAXIS  0,0. 005,0, XI 
70  YAXIS  0,0. 1,0, R 

80  LABEL  (*,1.5, 1.7, 0,1) 

81  PLOT  X1/2,R/2,1 

82  IF  A#1  THEN  85 

83  LABEL  (*)  "FRONT  FACE" 

84  GOTO  90 

85  IF  A# 2 THEN  88 

86  LABEL  (*)  "SIDE  FACE" 

87  GOTO  90 

88  LABEL  (*)  "BACK  FACE" 

90  FOR  Y=0.5  TO  R STEP  0.5 
100  PLOT  0,Y,1 

110  CPLOT  -5,-0. 3 
120  LABEL  (170)Y 
130  NEXT  Y 
140  PLOT  -7,R/2,1 
150  CPLOT  -7,-0. 3 

160  LABEL  (*,1.5, 1.7, PI/2,1)  "PRESSURE,  PSI" 

170  FORMAT  F4.1 

180  LABEL  (*,1.5, 1.7, 0,1) 

190  FOR  X=0.01  TO  Xl  STEP  0.01 
200  PLOT  X,0,1 
210  CPLOT  -3,-2 
220  LABEL  (240) X 

230  NEXT  X 

231  PLOT  Xl/2,-4,1 

232  CPLOT  -5,-0. 3 

233  LABEL  (*)  "TIME,  SECS" 

240  FORMAT  F6.2 

250  0=14 . 7*5/14* (P0/14 . 7) #2/ ( 1+1/7*P0/14 .7) 

260  IF  A#1  THEN  310 
270  X2=3*S/U 

280  Y2=Q*(1-X2/X1)  2*EXP (-2*X2/X1 ) +P0* ( 1-X2/X1 ) *EXP (-X2/X1 ) 

290  L1=X1 

300  GOTO  400 
320  X2=L/U 

330  Y2=-0.4*Q*(1-X2/X1/2)42  EXP(-2*X2/X3  ^2) +P0* ( 1- X2/X1/2 ) *EXP (-X2/X1/2) 

340  Ll=Xl+X2/2 
350  X3=(L+4*S)/U 

360  X2=L/U 

370  Y2=-0,4*Q*(1-(X3-X2)  A1)42*EXP(-2*(X3-X2)/X1) 

380  Y2=Y2+P0* (1- (X3-X2) /XI ) *EXP (- (X3-X2 ) /XI ) 


390  L1=X1+X2 

400  D=X1/D1 

410  PRINT  "TIME",  "PRESSURE" 

420  FOR  X=0  TO  LI  STEP  D 
430  GOSUB  A OF  470,540,620 
440  NEXT  X 
450  PEN 
460  STOP 

470  y=0.4*Q* (1-X/Xl)f 2*EXP(-2*X/X1)+P0* (1-X/Xl) *EXP (-X/X1) 

480  IF  X2  D THEN  510 
490  IF  X >=  X2  THEN  510 
500  Y=R* (l-X/X2)+y2*X/X2 
510  PLOT  X,y 
520  PRINT  X,Y 
530  RETURN 

540  N=0,4*Q* (l-(X-X2/2)/Xl)  2*EXP (-2* (X-X2/2) /XI) 

550  Y=N+P0* (1- (X-X2/2)/Xl) *EXP (- (X-X2/2) /XI) 

560  IF  X2  D THEN  590 
570  IF  X>=  X2  THEN  590 
580  Y=Y2*X/X2 

590  PLOT  X,y 
600  PRINT  X,Y 
610  RETURN 

620  Y=0. 4*Q* (1- (X-X2) /XI)  2*EXP (-2* (X-X2) /XI) +P0* (1- (X-X2) /XI) *EXP (- (X-X2) 

/XI) 

630  IF  X3  D THEN  680 
640  IF  X>=  X3  THEN  680 
650  y=0 

660  IF  X<X2  THEN  680 
670  Y=(X-X2)/(X3-X2)*Y2 

680  PLOT  X,Y 
690  PRINT  X,Y 
700  RETURN 

9000  DATA  1,0.05,40,1,0.5,2 
9990  END 


14 


TIME 

0 

1.25000E-03 

2.50000E-03 

3.75000E-03 

5.00000E-03 

6.25000E-03 

7.50000E-03 

8.75000E-03 

0.01 

0.01125 

0.0125 

0.01375 

0.015 

0.01625 

0.0175 

0.01875 

0.02 

0.02125 

0.0225 

0.02375 

0.025 

0.02625 

0.0275 

0.02875 

0.03 

0.03125 

0.0325 

0.03375 

0.035 

0.03625 

0.0375 

0.03875 

0.04 

0.04125 

0.0425 

0.04375 

0.045 

0.04625 

0.0475 

0.04875 

0.05 


PRESSURE 

2.057747834 

1.76201578 

1.466283727 

1.170551673 

0.874819619 

0.777923667 

0.736753283 

0.697168292 

0.659113619 

0.622536070 

0.587384259 

0.553608550 

0.521160990 

0.489995253 

0.460066578 

0.431331720 

0.403748894 

0.377277723 

0.351879192 

0.327515598 

0.304150507 

0.281748709 

0.260276177 

0.239700027 

0.219988476 

0.201110810 

0.183037342 

0.165739381 

0.149189198 

0.133359992 

0.118225860 

0.103761768 

0.089943520 

0.076747731 

0.064151801 

0.052133885 

0.040672875 

0.029748369 

0,019340650 

9.43066E-03 

0 
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APPENDIX  C 


Computer  Program  for  Generation  of  NASTRAN-DAREA 
Cards  for  Loaded  Grid  Points 

A FORTRAN  program  (LOAD) , for  use  on  the  Burroughs  B5700  was  devised  to 
generate  the  NASTRAN  DAREA  cards  for  a uniform  time-dependent  pressure  on 
the  front  or  back  surface  of  the  reflector.  The  program  reads  as  input 
the  GRID  cards  and  connection  member  cards  for  trianglular  or  quadrangular 
plates  on  the  surface.  These  are  preceded  by  a data  card  with  three  integers 
on  it,  separated  by  commas.  The  first  integer  is  a flag  which  equals  one 
(1)  for  a load  on  the  front  face.  The  second  integer  indicates  the  number 
of  grid  points  to  be  read,  while  the  third  indicates  the  number  of  plates. 

The  GRID  cards  follow  in  any  order  and  the  connecting  member  cards  follow 
the  grid  cards  in  any  order.  A listing  of  the  program  appears  on  the 
following  pages. 


i 


1 ONSITE 

6 FILE  13=HERB,  UNIT=READER 

10  REAL  RI,R(200) ,TI,T(200) ,ZI,Z(200) ,K,HI,KI,H4,AA,AB,AC,AD,L 

20  REAL  CA,CB,CD,CC,LR(200) ,LZ (200) 

21  REAL  LE 

30  INTEGER  NG ,GR,G (200) , NPL, A, B,C,D, A1 , B1 ,C1 ,D1 

31  INTEGER  LARRY (8) , BARRY (8) , HARRY (8) , GARRY (8) 

35  READ  (13,/)PJ,NG,NPL 

36  PRINT  51,PJ,NG,NPL 
40  DO  600  1=1, NG 

50  READ(13,800)LARRY,RI,TI,ZI 

51  51  FORMAT(3I3) 

52  PRINT  800, LARRY, RI,TI,ZI 
55  CALL  HAROLD  (LARRY, GR) 

60  G(I)=GR 

70  R(G(I))=RI 
80  T(G(I) )=TI*3. 141593/180 
90  600  Z(G(I))=ZI 

100  800  FORMAT  (8X , 8A1 , 8X , 3F8 , 3) 

110  801  FORMAT  (24X,4(8A1)) 

120  802  FORMAT  ("CHECK  GEOMETRY") 

130  DO  604  1=1, NPL 

131  READ (13, 801) LARRY, BARRY, HARRY, GARRY 

132  PRINT  801, LARRY, BARRY, HARRY, GARRY 

141  CALL  HAROLD  (LARRY, Al) 

142  CALL  HAROLD  (BARRY, Bl) 

143  CALL  HAROLD  (HARRY, Cl) 

144  CALL  HAROLD  (GARRY, Dl) 

150  IF  (Dl.NF.O)  GO  TO  603 
160  D1=A1 

170  603  IF  (R(A1) .EQ.R(Dl) ) GO  TO  606 

180  IF  (R(C1) .EQ.R(Bl)  ) GO  TO  602 
190  PRINT  802 
200  STOP 

210  602  IF  (T(A1) .NE.T(Bl) ) GO  TO  605 

220  A=B1 

230  B=A1 

240  C=D1 

250  D=C1 

260  GO  TO  601 

270  605  A=C1 

280  B=D1 

290  C=A1 

300  D=B1 

310  GO  TO  601 

320  606  IF  (T(A1) ,NE.T(B1) ) GO  TO  607 

330  A=A1  ■ 
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340  B=B1 
350  C=C1 
360  D=D1 
370  GO  TO  601 
380  607  A=D1 

390  B=C1 
400  C=B1 
410  D=A1 

420  601  K=(R(B)-R(A))*COS( (T(D)-T(A))/2) 

430  HI=(K**2+(Z(B)-Z(A))**2)**.5 

440  K1=(R(C)-R(D) )*COS( (T(D)-T(A) )/2) 

450  H4= (KI**2+(Z(C)-Z(A) )**2)**.5 

460  AA=R(A)*SIN((T(D)-T(A))/2)*HI 

470  AD=R(A) *SIN( (T (D) -T (A) ) /2) *H4 
480  AB=R(C) *SIN( (T(C)-T(A) )/2)*HI 

490  L= (R(B)-R(C) )* (R(A)*SIN( (T (D) -T (A) ) /2) / (R (B) -R (A) ) ) 

500  AC-HI* (R(C)*SIN( (T(C)-T(A) )/2)-L) 

501  AA-ABS (AA) 

502  AB=ABS (AB) 

503  AC=ABS(AC) 

504  AD=ABS(AD) 

510  CA= (AA+AB+AD)/6 
520  CB=  (AA+AB+AO/6 
530  CC= (AB+AC+AD)/6 
540  CD= (AC+AD+AA)/6 

546  K=ABS(K) 

547  HI=ABS(HI) 

550  608  LR(A)=LR(A)+CA+K/HI 

560  LR(B)=LR(B)+CB*K/HI 
570  LR(C)=LR(C)+CC*K/HI 
580  LR(D)=LR(D)+CD*K/HI 
590  LZ(A)=LZ(A)+CA*LE 
600  LZ(B)=LZ(B)+CB*LE 
610  LZ(C)=LZ(C)+CC*LE 
620  LZ(D)=LZ(D)+CD*LE 
630  604  CONTINUE 

650  804  FORMAT  ("DAREA  " , 2I8,.3X,  " 1 " ,4X,F8. 3 . 18, 3X , " 3"  , 4X  ,F8 . 3) 

670  DO  609  1=1, NG 

671  IF  (PJ.EQ.2)  GO  TO  610 

672  LR(G(I) )=-LR(G(I) ) 

673  GO  TO  609 

674  610  LZ(G(I) )=-LZ(G(I) ) 

675  PRINT  804,PJ,G(I) ,LZ(G(I) ) ,G(I) ,LR(G(I) ) 

680  609  PUNCH  804 , PJ ,G ( I) , LZ (G ( I ) ) ,G ( I ) , LR (G ( I ) ) 

690  STOP 

999  END 

2000  SUBROUTINE  HR (NAME, Q) 
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2010  DIMENSION  Q(8) 

2020  OS="00000" 

2030  DO  1 1=1,8 

2040  0T=CONCAT(0, 0(1) ,42,12,6) 

2050  IF(QT.E0.QS)GO  TO  1 
2070  NAME=NAME*10+OT 
2080  1 CONTINUE 

2090  RETURN 
2100  END 

10090  SUBROUTINE  HARDOL (LETTER, IV AL) 

10100  INTEGER  LETTER (8) , NUMS ( 10 ), BLANK 

10110  DATA  NUMS/lHO , IHl , 1H2 , 1H3 , 1H4 , 1H5 , 1H6 , 1H7 , 1H8 , 1H9/, BLANK/ IH/ 
10120  IVAL=0 
10122  DO  10  1=1,8 

10124  IF (LETTER (I) .NE. BLANK)  GO  TO  1000 
10126  10  CONTINUE 

10128  RETURN 

10130  1000  IF (LETTER (8) .NE. BLANK) GO  TO  2500 

10140  LSAVE= LETTER (8) 

10150  D01800  1=1,7 

10160  ia00LETTER(9-I)=LETTER(8-I) 

10170  LETTER (1)=SAVE 
10180  GO  TO  1000 
10190  2500  IDEC=0 

10200  2600  IF ( LETTER ( 8-IDEC) .LT. NUMS (1) ) GO  TO  3500 

10205  IF (LETTER (8-IDEC) .GT. NUMS (10) )GO  TO  3500 
10210  DO  100  1=1,10 

10220  IF(LETTER(8-IDEC).EQ.NUMS(I))GO  TO  3000 
10230  lOOCONTINUE 

10240  3000  IVAL=IVAL+ (I-1)*10**IDEC 

10250  IDEC=IDEC+1 
10260  GO  TO  2600 

10270  3500  IF(LETTER(8-IDEC) .EQ.BLANK)RETURN 

10280  PUNCH  200 

10290  200  FORMAT (17H0NVAL ID  CHARACTER) 

10300  RETURN 
10310  END 
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